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Liquid metals exhibit remarkable mechanical properties, in particular large surface tension and 
low viscosity. However, these properties are greatly affected by oxidation when exposed to air. We 
measure the viscosity, surface tension, and contact angle of gallium (Ga) and a eutectic gallium- 
indium alloy (eGaln) while controlling such oxidation by surrounding the metal with an acid bath 
of variable concentration. Rheometry measurements reveal a yield stress directly attributable to 
an oxide skin that obscures the intrinsic behavior of the liquid metals. We demonstrate how the 
intrinsic viscosity can be obtained with precision through a scaling technique that collapses low- 
and high-Reynolds number data. Measuring surface tension with a pendant drop method, we show 
that the oxide skin generates a surface stress that mimics surface tension and develop a simple 
model to relate this to the yield stress obtained from rheometry. We find that yield stress, surface 
tension, and contact angle all transition from solid-like to liquid behavior at the same critical acid 
concentration, thereby quantitatively confirming that the wettability of these liquid metals is due 
to the oxide skin. 



I. INTRODUCTION 

Liquid metals combine the largest surface tension of any fluid with high density and conductivity 1-3 . Several of 
these metals, in particular liquid gallium and its alloys, arc desired also for their low melting points and their wetting 
properties and have been used for thin film coatings 4 ' 5 , as moldable electronic interconnects 6 , or as penetrating 
contrast agents for the visualization of grain boundaries or micro-cracks 7-9 . However, the dynamical properties of 
liquid metals are significantly more complex than those of simple Newtonian liquids. At small applied stresses liquid 
gallium and its alloys show a solid-like elastic response 10,11 , and high-temperature molten metals under steady state 



shear have been found to exhibit non-Newtonian behavior 12 . It has been suggested 10 ' 11 that the solid-like properties 
are due to oxidation when the metal is exposed to air. It also has long been known 10 that gallium, which is highly 
wetting when oxidized, becomes non-wetting on glass when oxidation is prevented. A study with quantitative control 
of the oxidation conditions has not been performed so far. In particular, it has remained unclear if the desirable 
intrinsic properties of the pure liquid such as viscosity can be obtained along with wetting which is associated with 
oxidation. 

To address this issue, we present a systematic investigation of two liquid metals, gallium (Ga) and a eutectic 
gallium- indium alloy (cGaln). This study combines detailed measurements of the dynamic response, obtained by 
rheometry, with measurements of surface tension and contact angle. We control oxidation by surrounding the liquid 
metal in a bath of variable acid concentration. The acid bath not only prevents oxidation by mitigating contact with 
the air, but initiates a reduction process that can remove an oxide surface layer. By varying the acid concentration, 
we can precisely and reproducibly tune the balance between oxidation and reduction and map out the transition 
from solid-like to liquid behavior. This allows us to quantitatively analyze the role played by the oxide layer in 
determining the overall mechanical properties. 

The remainder of the paper is organized as follows. In Sec. II. we describe the materials and rhcometcr setup and 
measurement procedures, In Sec. III. A.l we present measurements of shear stress over time to demonstrate how 
oxidation and reduction at different acid concentrations result in time-dependence of shear stress at a fixed shear 
rate. In Sec. III. A. 2&3 we present steady state measurements of shear stress vs. shear rate, then introduce a new 
scaling technique to obtain the intrinsic liquid viscosity in the presence of non-Newtonian effects from rheometer data 
at higher Reynolds numbers. In Sec. III. B we show measurements of the effective surface tension using the pendant 
drop method and develop a model to identify the contribution of yield stress to drop size. Contact angle data are 
provided in Sec. III. C. This section also compares the yield stress, surface tension, and contact angle measurements 
as a function of acid concentration, demonstrating that the transition from solid-like to liquid behavior in all cases 
occurs at the same critical acid concentration. Sec. IV contains a summary and conclusions. 

II. EXPERIMENTAL 
A. Materials 

We purchased Ga (purity 99.99%) and Indium (In, purity 99.995%) for our experiment from Sigma-Aldrich, with 
the melting points around 30°C and 156°C respectively. We produced eGaln by mingling 77% Ga with 23% In. The 
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FIG. 1: (a) A picture of the rheometer set up. (b). A schematic illustration with dimensions. 

mixture was kept at a temperature of 200°C in a quartz tube for about two hours, while argon gas was cycled through 
the tube to inhibit oxidation. The liquid form of eGaln (melting point ~ 15°C) was obtained after cooling down 
the mixture to room temperature. All measurements of mechanical properties were performed at 50° C in the liquid 
state. At this temperature, Ga has a density of p — 6.095g/ml and nominal viscosity of r/i = 1.90 x 10 _3 mAT ■ s[l], 
and eGaln has a density of p — 6.25g/ml and nominal viscosity of r\i — 1.86 x 10~ 3 miV ■ s[10]. For the acid bath, we 
used hydrochloric acid (HC1) of varying concentration. 



B. Methods 

To measure the viscosity of liquid metals, we used an Anton Paar rheometer with a parallel-plate geometry. Fig.l 
shows the experimental setup, with plate diameter D = 25.13 mm and gap size d which was varied from 0.2 to 2 mm 
for different measurements. The rheometer measures the torque T exerted to rotate the top plate and the angular 
rotation rate uj. These are used to calculate the average shear stress at the plate edge r = 16T/irD 3 and average 
shear rate at the edge 7 = Duj/2d. The stress resolution is 0.03 Pa. The effective dynamic viscosity is r/ e // = r/j in 
the steady state. To contain the acid bath, a cup was constructed consisting of a clear cast acrylic cylinder with a 
31.75 mm inside diameter glued to a titanium plate of 0.51mm thickness that formed the bottom. The cup also helps 
contain the liquid metal, which due to their high surface tension with the tool surfaces do not stick. For this reason, 
we fill the liquid metal to a height larger than the gap size d so there is no space for the liquid to spill outward. 
Additionally, a 13.10 mm thick acrylic cylinder is attached to the tool to protect it from both the acid and gallium. 
At shear rates above ~ 10 4 s _1 , the surface shape of the liquid starts to change from a variety of instabilities, so we 
only report measurements at lower shear rates where the geometry of the liquid region remains fixed. 
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III. RESULTS 
A. Viscosity 

1. Oxidation and reduction on the surface 
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FIG. 2: Oxidation of liquid metals exposed to air in rheometer measurements. The time evolution of shear stress r for both 
Ga (blue) and eGaln (red) at j = 100s -1 (open symbols) and 7 = 1200s -1 (closed symbols). The increase of r with time at 
7 = 1200s -1 implies ongoing solidification by oxidation at the surface. Inset: plot of dr/dt vs. time. 

We first demonstrate the effect of the oxidation process on the rheology of liquid metals. In these experiments we 
measured the time dependence of shear stress r at constant shear rate 7 as oxidation and reduction processes may 
be taking place, under different initial conditions, both in air and in the acid bath, and for both Ga and eGaln. 

Initially, 2ml of oxidized liquid metal was directly placed on the bottom plate of the rheometer with a pipette and 
exposed to the air, and then rested for around 30s to allow the shape of the surface to fully relax. Such an oxidized 
droplet appears dirty and wrinkled on the surface, while the bulk appears pure if we slice open the skin. Care was 
taken to avoid any agitation or shear that could break the layer before the experiment. We then lowered the tool to 
a gap size d = 2 mm and started to shear the sample with a constant shear rate. Fig. 2 shows the time-evolution of 
shear stress for experiments at different shear rates. The existence of the solid oxide layer is seen to initially cause 
a large stress r ~ 10 2 Pa, which is two orders of magnitude larger than the expected viscous contribution of the 
bulk 1 ' 10 for 7 = 100 s . At this shear rate, r displays very slow increase with time, indicating that the sample 
is almost in equilibrium. However, at 7 = 1200s -1 , r increases linearly with time. It can be seen that the oxide 
skin is broken at the higher shear rate and the surface layer turns over so that the bulk material becomes exposed 
to the air and has a chance to react with the oxygen, creating more solid oxidized material to resist the rotation of 
the rheometer plate. This result can be explained by continuing solidification by oxidation happening on the surface 
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under shear. 

To demonstrate the role of the acid bath in reducing the oxidized surface layer, we performed a similar experiment 
in acid baths of varying concentrations of HC1. In these experiments, 2ml of previously oxidized Ga (or eGaln) was 
placed with a pipette into the acrylic cup (Fig.l) filled with HC1 solution. A constant shear rate was imposed at 
t = immediately after the liquid metal was placed in the acid. Fig. 3 (a) shows the time-evolution of the shear stress 
for 7 = 60s -1 for several different acid baths with acid concentrations c H ci- In contrast with the measurements with 
exposure to air, now the shear stress starts at low values, decreases with time, and eventually reaches a steady state. 
This suggests the acid bath not only inhibits further oxidation but actually reduces the pre-existing oxide layer, and 
eventually an equilibrium value is reached depending on concentration where oxidation and reduction processes are 
in balance. The equilibrium stress is seen to decrease with chci an d level off for chci ~ 0.5 M. Equilibrated liquid Ga 
in acid baths with higher chci appears to be very shiny, comparable to a mirror finish, which is expected for the pure 
liquid. The remaining stress in equilibrium can be attributed to hydrodynamic stresses in the bulk. This suggests 
that a pure liquid state of Ga can be reached by vigorous washing of the liquid metal in HC1 at a concentration 
chci ~ 0.5 M. 

The time-evolution of stress for additional measurements at different shear rates with chci = 0.5 M are shown in 
Fig. 3 (b). The rate of stress decrease is seen to increase with 7. This conforms what we observed for the oxidation 
measurements: faster shear increases the exposure of the bulk liquid metal to the surrounding environment and 
accelerates the corresponding reactions. To quantify the relationship between reduction rate and shear speed, we 
calculate (—d,T/dt) m as the slope of the steepest portion along each curve, and plot this slope vs Reynolds number 
Re = pd 2r y/i]i in Fig. 3 (c). There appear to be two linear scaling regimes, with a steeper slope at higher Re. We 
fit the two regimes by straight lines, whose intersection defines a transitional Reynolds number ReJ <~ 21 ± 2 above 
which the rate of stress increase by oxidation increases rapidly. We attempted a similar measurement for the shear 
rate dependence of the oxidation in air (Fig. 2). However, the oxidation in air happens much slower than reduction 
in the acid bath (by a factor of <~ 10 2 ) so that we could not resolve a similar transition for oxidation. In Sec. III. 
A. 3, we will show that inertial flow is seen for Re > Re' c , suggesting the enhanced reaction rate is the result of 
enhanced mixing from vortical flows at high Re. 



5 



(a) 



10 2 



:22SS°°°oooo , 



CO 

Q- 10 1 



10°, 



1.0M HCI 
0.5M HCI 
0.2M HCI 
0.1M HCI 
Reduction Process 



19Pa 
Equlibrium State 



v — m.-i "■■ Equlibrium State 

V V V A V V 1 f, Or, 



(b) 



Q_ 10 1 



10", 



200 400 

Time t (sec) 

1.2 



C HCI =0.5/W 

***•!;««:::"■—■■■——.■.— .................1 




200 400 

Time t (sec) 




FIG. 3: Reduction of the oxide layer for gallium in rheometer measurements with an acid bath, (a) Time evolution of shear 
stress r at -y = 60 s _1 for different acid concentrations, (b) Time evolution of shear stress r at different shear rates for 
C HC1 = 0-5 M. (c). Maximum rate of change of stress — (4?)m from panel (b) plotted against Reynolds number. Here, dashed 
lines represent two different linear scaling regimes, whose intersections defines a transitional Reynolds number, Re' c ~ 21 ± 2, 
above which the reduction process speeds up dramatically. 




FIG. 4: Shear stress r vs shear rate 7 in the steady state for (a) Ga and (b) eGaln. HCI bath concentrations are shown in the 
keys. Solid circles: liquid metals are exposed to air without an acid bath. Open symbols: Gallium in an Argon atmosphere 
without an acid bath, in which the gallium was first brought to equilibrium either in air (open blue circles) or a 1 M HCI bath 
(open yellow diamonds). Solid line: r oc -y 3 / 2 corresponding to intermediate Reynolds number flow. 
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2. Viscosity curves 



We next characterize the relation between shear stress r and shear rate 7 for Ga and eGaln at different atmospheric 
and bath conditions in the steady state. To achieve chemical equilibrium, the samples were pre-sheared before each 
measurement at 7 = 60s" 1 for around 600s according to Fig. 3. The steady state was confirmed by slowly ramping 
the shear rate first up then down, and checking that there was no hysteresis in the measured stress, except when the 
liquid metal is exposed to air where the oxidation is not reversible by shear. Figure 4 shows t vs 7 at for different 
acid bath concentrations chci at a gap size d = 2 mm. Data are shown for Ga in panel (a) and for eGaln in panel 
(b). When exposed to air, the liquid metals are seen to exhibit a yield stress, i.e, the shear stress goes to a non-zero 
value in the limit of zero shear rate. This indicates solid behavior from the oxidized layer, since at stresses below the 
yield stress, the material does not flow. This value is consistent with oscillatory rheometer measurements 11 . 

Adding the acid bath and increasing the HCL concentration reduces this yield stress. The yield stress plateaus at 
a minimum value for HC1 concentrations greater than about 0.5M for both metals. Along with the data of Figs. 2 
and 3, this confirms that an HC1 concentration greater than about 0.5M is enough to minimize the contribution of 
the oxidation to the shear stress. 

Another major feature of Fig. 4 is that the stresses collapse onto a single curve at high shear rates. This curve 
extends down to lower shear rates for higher acid concentrations where the yield stress is lower. This behavior 
suggests the total stress comes from the sum of two components: the yield stress which is dependent on the acid 
concentration, and a hydrodynamic contribution dependent on the shear rate. The latter component scales as 
t cx 7 3 / 2 , which differs from the laminar regime of a Newtonian fluid. We will show in the next section that this 
behavior is the hydrodynamic shear stress for intermediate Reynolds number flows where inertial effects play a role. 

So far we have used the acid bath to prevent oxidation, but an alternative method would be to use an inert 
atmosphere. However, Ref.l points out that argon gas environment can not protect the the surface of gallium from 
oxidation. To test this statement, we sealed a plastic bag around the rheometer and filled it with argon gas to 
replace the air. Argon gas was pumped into a larger cylinder surrounding the sample, and air was allowed to escape 
through a hole at the top. Liquid Ga samples were first washed in an acid bath to remove all oxidation, then placed 
within the rheometer in the argon atmosphere. A plot of r vs 7 for this sample is shown as open orange diamonds 
in Fig. 4(a). This sample behaves similar to those submerged in high-concentration acid baths. This test lasted well 
over lhr, over which time oxidation would have occurred to produce a significant yield stress if this experiment had 
been done in air. Therefore our results confirms that an inert atmosphere can prevent oxidation just as well as an 
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FIG. 5: Reynolds number collapse of stress measurements, (a). Viscosity rj e ff = r/j vs. shear rate 7 for liquid Ga in a 1.0M 
HC1 bath where there is no oxide layer. The gap size d is shown in the key (b). The rescaled shear stress T/{ryy) plotted 
against Re = pijd 2 /r/i for Ga, eGaln, and water. The data collapse in Re identifies two scaling regimes: the laminar regime 
at low Re where r oc 7 and an intermediate Re regime where r oc j 3 ^ 2 . The intersection of these two regimes identifies a 
critical Reynolds number Re c — 22 ± 3 for our measuring geometry. 

acid bath, which is in disagreement with claims in Ref.l. If instead a sample of liquid Ga that is already oxidized 
from exposure to air is placed in the rheometer with the argon atmosphere, the plot of r vs 7 shown as open blue 
circles in Fig. 4 (a) exhibits a large yield stress, similar to the case of a sample measured in air. This illustrates that, 
while an inert atmosphere can prevent further oxidation, it cannot reduce a pre-existing oxidation layer like an acid 
bath can. 

3. Reynolds number scaling 

In this section, we demonstrate that rescaling r(j) using the Reynolds number can be used to precisely measure 
the low intrinsic viscosities of liquid metals. The effective viscosity of a fluid is defined as ?7 e // = t/'J. When there 
is an oxide layer, 7y e // is dominated by the yield stress at low shear rates. To avoid this effect, we concentrate 
on measurements in a 1.0M HC1 bath, and plot r) e ft vs 7 for Ga in Fig. 5 (a) for different gap sizes d. In the 
limit of low shear rates, for the smallest gap sizes, ?7e// is constant, corresponding to the intrinsic bulk viscosity 
T]i ~ (1.82 ±0.21) x 1CT 3 Pa-s. 

This Newtonian regime is difficult to resolve in a rheometer for low viscosity fluids because the stress resolution 
limit leads to relatively large uncertainties at low shear rates, and many rheometers would not be able to resolve this 
regime. To our knowledge this is only the second direct rheometer measurement of the intrinsic viscosity of a liquid 
metal 12 . The Newtonian regime may also be unresolvable in the presence of non-Newtonian effects such as a yield 
stress; we have already omitted our data with a yield stress from Fig. 5, but a similar scaling technique could also 
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be applied to such data. In fact, we were only able to resolve the Newtonian regime at small gap sizes. At higher 
shear rates or larger gap sizes, the Reynolds number is larger so inertial effects increase energy dissipation and r) e ff 
increases with 7 as seen in Fig. 5 (a). 

To confirm that the observed increase in rj e f f with shear rate is an inertial effect, we rescale the data in Fig. 5 (a) 
non-dimensionally by dividing the stress by the viscous contribution and r/rnj and converting the shear rate to the 
Reynolds number Re = pcPj/r]i and replot it in Fig. 5 (b). It is now seen that the data for different gap sizes collapse 
onto each other, confirming that this dimensionless scaling is appropriate to describe the fluid flow. Particularly, at 
small gap sizes (light blue and pink points in Fig. 5), the plots extend along the constant viscosity region with a large 
range of shear rate, which thus can be used to determine viscous regime. By contrast, data at large gap sizes (red, 
dark blue and green points) is at a higher Reynolds number and thus larger stresses and a better stress resolution. 
Since this rescaling depends only on Newtonian fluid properties, it should apply to any Newtonian fluid in the same 
flow geometry. To test this, we also show measurements for water and eGaln in the same flow geometry in Fig. 5 
(b). The data for each liquid collapses onto the same curve, confirming that this rescaling applies to all Newtonian 
fluids in the same flow geometry. 

Since rheometer measurements have a better resolution in the high Reynolds number regime, we can take advantage 
of the universal scaling collapse to infer the intrinsic liquid viscosity from this regime by calibrating with a liquid of 
known viscosity. In the high Reynolds number regime (30 < Re < 6000), the dimensionless stress is found to follow 
the empirical relation t/(jjj7) oc Re 1 ! 2 for a fit parameter o, as shown in Fig. 5 (b), which is similar to other flows 
in a similar range of Re as a result of increasing inertial forces 13 . Combining this scaling regime with the viscous 
scaling r = 77/7 in the limit of low Re gives a relationship where the two scaling laws meet at a critical Reynolds 
number Re c defined by 




(1) 



Rearranging this gives a relationship gives an expression for T(j) 




(2) 
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Fitting this equation to r(j) gives the ratio between r]i/Re c . This process can be applied to a Newtonian fluid in 
the same geometry with known viscosity r\i to obtain Re ci so that rji can be determined for other fluids. 

Since the empirical inertial scaling is only valid in a limited range of Re, it is necessary that these fits are done in 
the same range of Re for different fluids. In our case, we measured water with known viscosity and simultaneously 
fit data for Re < 20 to the viscous relation r = rjij and the data for 30 < Re < 1000 to Eqn. 2 to obtain 
Re c — 22 ±3. The uncertainty is one standard deviation in the statistical uncertainty on the intersection point of the 
two simultaneous fits. Fitting the data for Ga and eGaln to the same functions in the same range of Re and using 
the known Re c gives f] Ga = (1.92 ±0.08) x 10~ 3 Pa- s for Ga and r/ cGaIn = (1.86 ±0.09) x 10" 3 Pa-s for eGaln, which 
are consistent with literature values 1,10 . The uncertainty here indicates one standard deviation in the statistical 
uncertainty on the viscosity from the simultaneous fit of the data to the two regimes. This result demonstrates that 
the intrinsic liquid viscosities can be measured with higher precision by using a Reynolds scaling collapse when stress 
resolution limits of rheometers make low-Re measurements difficult. 

This critical Reynolds number Re c = 22 ± 3 where inertial forces begin to contribute to energy dissipation during 
shear is consistent with the critical Reynolds number Re c ' ~ 21 ± 2 where the reduction rate of the oxide layer in 
acid bath increased (Sec. III. A. 1). This is not surprising, since as flows transition from laminar to non-laminar 
with increasing Re due to increasing inertial forces, vortices form which enhance transport and mixing of fluid. This 
mixing can expose more liquid metal to the surface to enhance oxidation or reduction rates at the same time. 



B. Surface Tension 



In this section, we show that the oxide layer on a liquid metal can affect the effective surface tension in addition to 
the effective viscosity. We performed pendant drop experiments, a standard technique for measuring surface tension. 
Droplets of liquid metal were formed on the end of a nozzle of diameter D n = 5.4 mm and 1.96 mm. Backlit images 
of Ga droplets of different sizes submerged in water are shown in Fig. 6 (a) . The effective surface tension j e f f can 
be obtained by fitting the equilibrium shape of the liquid metal droplet with the Laplace equation: 



AP + Apgz = leff {±- + ±-). (3) 

ill 

Here, R\ and i?2 are the two principle radii of curvature. Ap is the density difference of the two phases. The left- 
handed side of Eqn. 3 represents the pressure difference across the local interface, where AP is the pressure difference 
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FIG. 6: Effective surface tension measurements of Ga and eGaln using the pendant drop method, (a). Pendant droplets of 
liquid Ga submerged in water, (b). Measured y e g vs drop diameter I for Ga and eGaln in different HC1 baths concentrations 
shown in the key. Solid lines: fits of ^eff = tqI/A + 70, which give the pure liquid surface tension 70 (dashed lines) and the 
contribution from the oxide layer to, which is found to be equivalent to the yield stress r y from Fig. 4. 



at a selected reference plane and Apgz stands for the hydrodynamic pressure. Here, we fit Eqn. 3 numerically using 
the Axisymmetric Drop Shape Analysis (ADS A) technique 14-16 to obtain the effective surface tension j e ff of Ga 
and eGaln. A complete quantitative description of the ADS A method is given in the Appendix. We plot j e f f in 
Fig. 6 (b) for Ga and eGaln with different drop diameters I and different HC1 bath concentrations. For the water bath 
{chci = 0), the fit value of 7 c ff increases with I, which differs from a pure liquid where 7e// should be independent 
of I. In contrast, in the 1.0 M acid bath, the values of 7 e // shown in Fig. 6 (b) are seen to be independent of I, This 
suggests that the increase in 7 e // with drop size in water is due to the oxide layer. We also note that in the 1.0 M HC1 
bath, no stable drops were formed with I > 8 mm. In a pure liquid the largest drop size is limited by gravitational 
forces overcoming surface tension to pull the drop apart. Thus the oxide layer also allows larger droplets to remain 
stable. 

To understand the role of the oxide skin in the pendant drop measurements, we propose that the oxide skin 
provides an additional stress To along with the stress from surface tension to hold the drop together. For a simple 
approximation, we write a spatially averaged Laplace equation (Eqn. 3) with the effective surface tension split into 
pure liquid surface tension and yield stress terms 
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AP= 2 f. 2 | + , , (4) 



where 70 is the surface tension of the pure liquid and R is the average radius of curvature of the droplet. Estimating 
this as R rts 1/2 yields an expression for the effective surface tension as a function of drop size 



7off - ^ + 7o ■ (5) 

To test this size-dependence, we fit Eqn. 5 to the data in Fig. 6 (b), varying both r and 70 as fit parameters. For Ga 
in water, the slope gives t = 172 ± 12Pa, which is consistent with the yield stress measured in the rheometer without 
an acid bath r y = 189 ± 12Pa (Fig. 4). Similarly, for eGaln we obtain r = 119 ± 13Pa while r v = 108 ± 9Pa. These 
agreements indicate that the yield stress r y of the oxidized liquid metals is responsible for the larger, size-dependent 
effective surface tension in the pendant drop measurements. In the limit of I — or large chci, the effective surface 
tension in Eqn. 5 reduces to the surface tension of the pure liquid. Fits of Eqn. 3 yield 7 (Ga) = 695 ± 23mN/m 
and 70 (eGaln) = 445 ± 19mN/m, which arc in a excellent agreement with the pendant drop measurements in the 
acid bath shown in Fig. 6(b), as well as with literature values (Ga: 715mN/m; eGaln: 444mN/m ) 2 > 3 > 10 . These 
quantitative agreements confirm that we can rewrite Eqn. 5 as 7 e a = r y l/A + 70, where the effective surface tension 
in mechanical measurements is the linear sum of contributions from the liquid surface tension and the yield stress. 

We note that since the yield stress due to the oxide layer is a surface effect, it is not size-independent like an 
intrinsic property, but rather it should scale as r y = k/(l/2) — 2k /I where k is the effective stiffness of the oxide 
layer. 10 ' 11 This would change Eq. 5 to -f e ff ~ k/2 + -f . In our experiments, the gap size in the rheological test and 
the pendant drop size were similar, so it was not necessary to correct for the size dependence of the yield stress. 

C. Wetting 

We quantify the wetting properties of Ga and eGaln on a glass microscope slide using the contact angle 6 C . The 
pre-oxidized Ga and eGaln droplets were placed on the slide with a pipette, then allowed to equilibrate with the 
surrounding acid for 30 min. To avoid variations in the hydrostatic pressure from experiment to experiment, we 
used fixed volumes of liquid metal (250 y(il) and acid (500 ml) in each experiment. In each case, the droplets were 
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FIG. 7: Measurement of the contact angle 6 C for droplets of Ga on a glass surface submerged at different HC1 concentrations 
labeled in the images. 

pressed down against the glass surface using a spoon. For cjjci > 0.2M, the drops spontaneously bounced back to 
a spherical shape after removing the external force. However, at lower HC1 concentrations -where the oxide skin 
also contributes significantly to the yield stress and pendant droplet measurements - the droplet does not bounce 
back and the contact angle can vary significantly in repeated measurements. Figure 7 shows images of these Ga 
droplets on glass surfaces after this procedure in acid baths with different chci- We measure the contact angle 
between the droplet and glass surfaces on both sides of the image several times and take the average as 6 C . For each 
concentrations, and the standard deviation gives the uncertainty. These values of 9 C are plotted in Fig. 8 (e) and (f) 
for Ga and eGaln, respectively, in acid baths of different concentrations. The contact angle is seen to increase with 
HC1 concentration up to a critical concentration where 9 C reaches about 180°; this is the maximum possible contact 
angle, corresponding to a perfectly non-wetting droplet sitting on top of the glass surface. This suggest that pure 
liquid gallium is perfectly non-wettting on glass, and only wets the glass at lower acid concentrations where some 
oxide layer remains. 

We can compare the critical HC1 concentration found for the wetting on a glass surface to the effective viscosity 
and surface tension measurements. The yield stress r y from rheometer measurements (Fig. 4) is shown vs Chci in 
Fig. 8 (a) and (b) for Ga and eGaln, respectively. Similarly, the effective surface tension 7 e // is shown vs Chci in 
Fig. 8 (c) and (d) for Ga and eGaln, respectively. In each case, the scaling behavior of the physical parameter (r„, 
Jeff, or $c)j changes at a critical HC1 concentration of chci — 0.2 M for Ga and cjjci = 0.1 M for eGaln. Above 
this critical concentration, the material behaves as a pure liquid metal with no yield stress. The values for viscosity 
and surface tension arc in agreement with literature values, and the liquids arc perfectly non-wetting on glass. Below 
the critical HC1 concentration, a yield stress appears, the effective surface tension appears higher, and the materials 
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FIG. 8: Summary of physical parameters as a function of HC1 concentration chci for Ga (a,c,e) and eGaln (b,d,f). (a,b) 
Yield stress r y . (c,d) effective surface tension 7 e //. Different points correspond to different pendant drop sizes, where upper 
points generally correspond to larger pendant drops. (e,f) contact angle 8 C . Vertical dashed lines: critical HC1 concentration 
chci = 0.2M for Ga and 0.1M for eGaln which separates pure liquid behavior from oxide skin effects. 



are able to partially wet glass. The fact that these mechanical properties are all controlled by the same critical acid 
concentration confirms that they are related in oxidized liquid metals. In particular, this quantitatively shows that 
the wetting ability of Ga and eGaln is dependent on oxidation. 



IV. SUMMARY AND CONCLUSION 



Our results demonstrate that the mechanical properties of Ga and eGaln can be finely and reversibly controlled 
by adjusting the level of oxidation with an acid bath. Oxidation of liquid metals exposed to air leads to an increase 
in the shear stress (Fig. 2) and also produces a large yield stress (Fig. 4). This yield stress can be eliminated 
by surrounding the metal with acid, which prevents further oxidation from contact with air and reduces existing 
oxidation from the surface (Fig. 3). High shear rates enhance the oxidation and reduction process (Figs. 2, 3). We 
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showed that these chemical reactions accelerate once a critical Reynolds number is exceeded that signals the onset 
of enhanced mixing due to vortical flows (Fig. 5). 

The limited stress resolution of typical rheometers makes it difficult to resolve the viscosity of low- viscosity fluids 
in the laminar flow regime. To determine the intrinsic viscosity of liquid metals from rheometry data we introduced a 
new technique (Fig. 5). It takes advantage of the fact that, for fixed measuring geometry, all Newtonian fluids must 
exhibit the same scaling relation between dimensionless stress and Reynolds number. Therefore, a calibration with 
a fluid of known viscosity is sufficient to determine the viscosity of other Newtonian fluids in the same experimental 
set-up from data at higher Reynolds number. 

Pendant drop measurements showed that an oxide layer increases the surface tension compared to that of the pure 
liquid (Fig. 6). The resulting effective surface tension can be expressed as 7 e g w T y l/A + j , a linear combination of 
the pure liquid surface tension 70 and a contribution from the same yield stress t v from rheometry (Fig. 4). This 
oxide-dependent increase in effective surface tension may explain the significant variation in values of the surface 
tension for Ga reported in the literature 1 12 . 

To characterize wetting, we measured the contact angle of liquid metals on a glass surface, and found that the oxide 
layer enhances wetting by reducing the contact angle. Without the oxide layer, liquid Ga and eGaln are perfectly 
non- wetting (contact angle 180°) on a glass surface (Fig. 7). Remarkably, for effective surface tension, effective 
viscosity, and contact angle, the scaling behavior changed at the same critical concentration of the acid bath (Fig. 
8). Below this critical acid concentration the oxide layer resulted in a yield stress, large effective surface tension, 
and improved wetting, while above the critical concentration the intrinsic liquid properties were recovered. This 
provides a quantitative connection between oxidation and wetting properties of liquid metals. It also demonstrates 
a relationship between effective surface tension and viscosity that is unusual for Newtonian fluids. 
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FIG. 9: ADSA method, (a). The apex coordinate system used in ADSA. Here, the coordinates (x,z) are expressed as 
parametric form in term of s. Therefore, the local slope of (x(s), z(s)) is given by tan^ = x' (s) / y' (s) . (b). APn vs. cross 
diameter I for Ga (eGaln) in the solution with different acid concentrations, (c). 9o vs. cross diameter I for Ga (eGaln) in 
the solution with different acid concentrations. 



Appendix: Axisymmetric Drop Shape Analysis (ADSA) 



Here we summarize how to apply the ADSA method, as described in Ref. 17, to our sample. The geometry in 
our experiment is an axisymmetric pendant droplet, which is described by the apex coordinate system displayed in 
Fig. 9 (a). Mathematically, the shape of the axisymmetric drop can be presented by a parameter s with constant 
differential arc length dl = yj[x'{s) 2 + y'{s) 2 ] ds = Lds (0 < s < 1). The arc length between points O and A is then 
L = L ^/[x'{s) 2 + y'(s) 2 } ds. Using the curvature formula for parametric coordinates, the expression for Ri and i?2 
becomes 



1 z'{s) ' 2 x'{s)y"(s)-y'(s)x"{ S y [ > 

where the sign of i?2 is determined by the normal direction. For numerical consideration, the physical parameters 
in Laplace's Equation (Eqn. 3) can be scaled by nozzle radius R, surface tension 7 and density difference Ap. As a 



1G 



consequence, the Eq. 3 can be rewritten equivalently as coupled differential equations 



x"(s) = z"{s){2LP - LB z(s) + 



r(s) 



(7) 



z"{s) = -x"{s)(2LP - LB z(s) + 



r(s) 



(8) 



where x = x/R, z = z/R, L = L/R, P = AP R/(2j) and B Q = ApgR 2 /^f is the gravitational Bond number. In 
order to compute Eqs. 7 and 8 numerically, the boundary condition for contact angle 9o and pressure difference Ap 
at the edge of nozzle (point A in Fig. 9 (a)) need to be determined. Further, we utilized a Mathematica code of the 
ADS A method to search for the best 9 0} AP and surface tension 7, so that the fitted drop shape has the least mean 
square error. To test the validity of this method, we applied it to pendant droplets of water, where we obtained 
Iwater = 71.3 ± 2.3mN/m (literature value: 72.8 mN/m). 

Another benchmark can be carried out by checking the value of fitted boundary condition in different measurements. 
Fig. 9 (b) & (c) shows the size-dependence of AP and #o, respectively. For both Ga and eGaln, AP collapse onto 
constant values for different acid concentrations: AP (Ga) <~ 900Pa and AP (eGaln) ~ 700Pa. These results are 
consistent with the hydrostatic pressure in the syringe with liquid height « 13 mm. On the other hand, 9 does 
depend on I, especially when the drop size becomes large (I > 6mm). This transition may be determined by the 
capillary length where gravitational stresses overcome surface tension. Here, it is important that #o is not the same 
as 9 C measured in Sec. III. C. At the contact point, the liquid metal is actually touching the nozzle edge instead of 
a solid plane. Therefore, the liquid is pinned around the edge and can freely vary with the shape of the droplet. 



[1] K. E. Spells, "The determination of the viscosity of liquid gallium over an extended range of temperature" Proc. Phys. Soc. 
48, 299(1936) 

[2] S. C. Hardy, "The surface tension of liquid gallium", J. Cryst. Growth 71, 602-606 (1985) 

[3] B. B. Alchagirov, A. G. Mozgovoi, "The surface tension of molten gallium at high temperatures", High. Temp. 43, 
791-792(2005) 

[4] J. Ciosek, J. A. Dobrowolski, G. A. Clarke and G. Larframboise, "Design and Manufacture of All-Dielectric Nonpolarizing 

Beam Splitters", Appl. Optics. 38, 1244-1250(1999) 
[5] M. Ruckh, D. Schmid, M. Kaiser, R. SchfHer, T. Walter and H.W. Schock, "Influence of substrates on the electrical 

properties of Cu(In,Ga)Sc2 thin films", Sol. Energ. Mat. Sol. C.41-42, 335-343 (1996) 
[6] R. C. Chiechi, E. A. Weiss, M. D. Dickey and G. M. Whitesides, "Eutectic Gallium-Indium (EGaln): A Moldable Liquid 



17 



Metal for Electrical Characterization of Self- Assembled Monolayers", Angew. Chem. 120, 148-150 (2008) 
[7] W. Ludwig, D. Bellet, "Penetration of liquid gallium into the grain boundaries of aluminium: a synchrotron radiation 

microtomographic investigation", Mat. Sci. Eng. A-Struct. 281, 198-203 (2000) 
[8] R.C. Hugo, R.G. Hoagland, "The kinetics of gallium penetration into aluminum grain boundaries - In situ TEM obser- 
vations and atomistic models", Acta. Mater. 48, 1949-1957 (2000) 
[9] W. Ludwig, J.Y. Buffire, S. Savelli, P. Cloetens, "Study of the interaction of a short fatigue crack with grain boundaries 
in a cast Al alloy using X-ray microtomography" , Acta. Mater. 51, 585-598 (2003) 
[10] M. D. Dicky, R. C. Chiechi, R. J. Larsen, E. A. Weiss, D. A. Weitz, G. M. Whitesides, "Eutectic Gallium-Indium (EGaln): 
A Liquid Metal Alloy for the Formation of Stable Structures in Microchannels at Room Temperature" , Adv. Funct. Mater. 
18, 1097-1104 (2008) 

[11] R. J. Larsen, M. D. Dickey, G. M. Whitesides, D. A. Weitz, "Viscoelastic properties of oxide-coated liquid metals", J. 
Rheol. 53(6) 1305-1326 (2009) 

[12] M. Jeyakumar, M. Hamed, S. Shankar, "Rheology of liquid metals and alloys", J. Non-Newtonian Fluid Mech., 166(2011), 
831-838 

[13] R. J. Donnelly, N. J. Simon, "An empirical torque relation for supercritical flow between rotating cylinders, " , J. Fluid 
Mech. 7(3), 401-418 (1960) 

[14] Y. Rotenberg, L. Boruvka, A. W. Neumann, "Determination of Surface Tension and Contact Angle from the Shapes of 

Axisymmetric Fluid Interfaces", J. Colloid. Intcrf. Sci. 93, 169 (1982) 
[15] S. H. Anastasiadis, J. -K. Chen, J. T. Koberstein, A. F. Siegel, J. E. Sohn, J. A. Emerson, "The Determination of 

Interfacial Tension by Video Image Processing of Pendant Fluid Drops", J. Colloid. Interf. Sci. 119, 55, (1986) 
[16] M. G. Cabezas, A. Bateni, J. M. Montanero, A. W. Neumann, "Determination of Surface Tension and Contact Angle 

from the Shapes of Axisymmetric Fluid Interfaces without Use of Apex Coordinates", Langmuir 22, 10053-10060 (2006) 
[17] J. C. Burton, F. M. Huisman, P. Alison, D. Rogerson, P. Taborek, "Experimental and Numerical Investigation of the 

Equilibrium Geometry of Liquid Lenses", Langmuir 26(19), 15316-15324 (2010) 



18 



